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We have investigated temperature-dependent electronic structures of Na 2 lr 03 to unravel its in¬ 
sulating nature. Employing the combined scheme of the density-functional theory (DFT) and the 
dynamical mean-field theory (DMFT), we have shown that the insulating state persists even above 
the Neel temperature (Tjv), which reveals that Na 2 lr 03 is classified into a Mott-type insulator. The 
measured photoemission spectrum in the paramagnetic (PM) state is well described by the electronic 
structure obtained from the DFT-I-DMFT for the insulating state above Tjv. The analysis of optical 
conductivity, however, suggests that the non-local correlation effect is also important in Na 2 lr 03 . 
Therefore, Na 2 lr 03 is not to be a standard Mott insulator in that the extended nature and the 
non-local correlation effect of Ir 5d electrons are important as well in describing its electronic and 
magnetic properties. 

PACS numbers: 75.47.Lx, 71.20.-b, 71.70.Ej 


Identifying the insulating nature of transition metal 
oxides has been a central and long-standing subject in 
modern condensed matter physics. [l| Recent attention 
has been paid to 5c? transition metal oxides, Sr2lr04 and 
Na2lr03, whether they belong to Mott-type or Slater- 
type insulators. The strong spin-orbit coupling (SOC) 
and the rather weak Coulomb correlation of 5c? electrons 
are known to be two essential ingredients in determin¬ 
ing the ground state physics of Sr2lr04 and Na2lr03.[l- 
2l| Despite intense studies on the role of interplay be¬ 
tween Coulomb interaction and SOC, however, there has 
been no consensus on the nature of their insulating states 
yet. For example, on the insulating nature of Sr2lr04, 
there exist two contradictory reports, Mott insulator[6j 
vs. Slater insulator. Q A marginal Mott insulating state 
was also proposed, in which the insulating state above 
the Neel temperature {Tn= 2 A 0 K) was attributed to 
the presence of short range antiferromagnetic (AFM) 
correlation. im For Na2lr03 too, which is a system 
of our present interest, there have been debates on its 
insulating nature. [lll - [l4 [I7l - l2lj| 


Na2lr03 exhibits insulating state at room tempera¬ 
ture (r), well above rAr=15 K. 0 [0 The paramag¬ 
netic (PM) state with Curie-Weiss susceptibility behav¬ 
ior was confirmed up to T = 500 K.[^ The zigzag-type 
AFM ordering occurs below Tv.(^.l2^ To explain the 
AFM insulating ground state of NaMrOs, both the Mott- 
type 0, 0 and Slater-type [0 [isl [0 I0| mechanisms 
of the metal-insulator transition were invoked. In the 
former, the good basis to describe the local electronic 
structure is derived by consideration of the on-site atomic 
SOC and the crystal held. This basis is usually very close 


24| 


to the relativistic Jeff= 1/2 orbital. S 0 0 12 
Superexchange-based formalism is used to describe the 
zigzag-type AFM ordering. 0 , 0 | The observed PM in¬ 
sulating state well above Tv supports the Mott-type 
mechanism. ll| In the Mott-type mechanism, however, 
very long range magnetic interaction or very high en¬ 


ergy excitation is required to describe the zigzag-type 
AFM ordering, which raises the question on the valid¬ 
ity of superexchange-based formalism. 0 On the other 
hand, in the Slater-type mechanism, the quasi-molecular 
orbital (QMO) has been suggested as a good basis with 
consideration of the strong anisotropic hybridization be¬ 
tween 5c? orbitals of neighboring Ir atoms. 0 ] The zigzag- 
type AFM ordering occurs as a consequence of the energy 
gain from the gap opening at the zone boundary. 0 211 
The Coulomb correlation of 5c? electrons is considered 
just to enhance the band gap. Larger extension of 5c? 
orbital and multi-peak features in photoemission spec¬ 
trum (PES) and optical conductivity seem to support 
the Slater-type mechanism. However, in the Slater-type 
mechanism, it is hard to explain the observed PM insulat¬ 
ing state well above T^. Therefore, the issue is whether 
the local AFM ordering is essential in describing the in¬ 
sulating state of Na2lr03 or not. 


In this letter, we have investigated T-dependent elec¬ 
tronic structures and magnetic properties of Na2lr03, 
using the combined scheme of the density functional 

theory (DFT) _ and the dynamical mean-held theory 

(DMFT),[0, 0| and unraveled its insulating nature. 
We have shown that Na2lr03 at room T has the elec¬ 
tronic structure of the PM insulating state, revealing 
that Na2lr03 is a Mott-type insulator. Highly incoher¬ 
ent PES for Na2lr03 in the PM state is well described by 
the DMFT incorporating the local dynamical correlation. 
However, we have also found evidences that Na2lr03 de¬ 
viates from standard Mott insulators due to extended na¬ 
ture of Ir 5c? orbitals. The onset of the zigzag-type AFM 
ordering induces the signihcant redistribution of charge 
and spin densities with respect to those of Jeff= 1 / 2 , or¬ 
bital. This feature implies that, even though Na2lr03 is 
classihed into a Mott-type insulator, the superexchange- 
based theory assuming the rigid charge degree of freedom 
needs to be refined to elucidate the magnetic ordering in 
Na2lr03. Moreover, optical conductivity, which corre- 
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FIG. 1: (Color online) Ir 5d partial density of states (PDOS) 
of Na 2 lr 03 from the DFT (a), and the DFT+DMFT at 
T = 300 K (b). Red solid, black dashed, and purple dotted 
lines represent jeff=ll2, je//=3/2, and eg projected PDOSs, 
respectively. See the text for numerical ye//=l/2, je//=3/2, 
and eg basis, (c) The spectral function from the DFT+DMFT 
at T = 300 K. High symmetry points of of Brillouin zone are 
presented in the Supplemental Material. [33] 


spends to two-particle property, is not well described by 
the DMFT incorporating only the local dynamical corre¬ 
lation. Thus, the additional non-local spatial electronic 
correlation is expected to be important to describe opti¬ 
cal conductivity in Na 2 lr 03 . 

We have performed the fully charge self-consistent 
DFT+DMFT calculations based on the projection¬ 
embedding scheme. [^ . Correlated Ir bd electrons 
were treated dynamically by the DMFT local self-energy, 
while s and p electrons were treated on the DFT level. 
The SOC in the Ir hd orbital is included in all the calcula¬ 
tions. For Coulomb interaction parameters, U and J, we 
employed U =3.5 eV and J=0.8 eV for the DFT+DMFT. 
These parameters yield the band gap and PES spectrum 
in good agreement with the experiment. The DFT part 
calculation was done by using the full-potential linearized 
augmented plane wave (FLAPW) band method. 3^, Slj 
For the DFT+DMFT calculation, T = 300 K was cho¬ 
sen, and the PM state is assumed. We adopted experi¬ 
mental crystal structure of Na 2 lr 03 .Q We used numer- 



FIG. 2: (Golor online) Ir 5d PES at T = 130 K (black 
dashedl [lll| is compared with calculated PDOSs in the 
DFT+f/ {U = 2.75 eV) and the DFT+DMFT (T = 300 K). 
Blue dotted and red solid lines represent the DFT+I/ and the 
DFT+DMFT PDOSs, respectively, which are broadened with 
a Gaussian function (0.10 eV FWHM). The Fermi level is ad¬ 
justed in the case of DFT+f/ to £t the experiment. Red solid 
and black dashed arrows point at two dominating peaks in 
the DFT+DMFT and the measured PES, respectively. Blue 
dotted arrows point at extra peaks in the DFT+f/, which be¬ 
come signihcantly broadened in the DFT+DMFT due to the 
correlation-induced incoherency. 


leal jeff=lf2, jejj=3/2, and bases for Ir 5d elec¬ 
trons in the DFT+DMFT, which diagonalize the infi¬ 
nite frequency hybridization function A(a; = oo) in the 
local impurity Green’s function of the DMFT result. 
As shown below, these jeff bases are consistent with 
relativistic Jpff p rojectors for t 2 g orbitals in the cubic 
ITjl In comparison, we have also performed 


If 


symmetry. 

the DFT and DFT+f/ calculations. 32| For the DFT+f/, 
we used U = 2.75 eV and J = 0.6 eV, following the 
existing calculation.!^ We assumed, for the magnetic 
structure, the observed zigzag-type AFM ordering. fsll^ 
See the Supplemental Material (SM) for computational 
details.!^ 

Figure [T])a) and (b) presents partial density of states 
(PDOS) of je//=l/2, je//=3/2, and eg states in the DFT 
and the DFT+DMFT, respectively. The PDOS in the 
DFT in Fig. 1(a) shows six peaks in the t 2 g part. Due 
to anisotropic solid environment in Na 2 lr 03 , the effec¬ 
tive orbital degeneracy in the t 2 g manifold is seen to be 
well lifted. Mazin et al. 13| interpreted these six peaks as 
the QMOs resulting from strong intra- and weak inter¬ 
hexagon hybridizations. On the other hand, in the two 
uppermost t 2 g-driven bands, jeff=1/2 orbital (red solid) 
has main contribution with a small amount of je //=3/2 
(black dot) contribution. These QMO-driven six peaks 
and dominant contribution from jeff=1/2 near Ep are 
consistent with previous experimental and theoretical 
reports. 1^ 17, The metallic PDOS in the DFT is not 
consistent with the insulating ground state of Na 2 lr 03 . 

The PDOS in the DFT+DMFT in Fig. [Ilb) shows 
clear Mott gap of ^^400 meV size in the jeff =1/2 states 


near Ep, as is consistent with experiment. [llj Note that 
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(a) Crystal Structure (b) DFT (NM) 



FIG. 3: (Color online) (a) Crystal structure of Na 2 lr 03 (side view in the left and top view in the right).Blue, yellow, and 
red balls correspond to Ir, Na, and O atoms, respectively. Numbers (1—6) in the top view are given to Ir atoms, for which 
spin and charge densities are plotted in (b), (d), (e)-(f). (b) The charge density plot of the unoccupied Ir t 2 g band of the 
nonmagnetic (NM) state in the DFT. Red balls represent oxygen atoms, (c) Projection (Z) between numerical je// = l/2 bases 
in the DFT and the DFT+DMFT at cj=0, and a;=oo with variation of U. (d) The spin density plot of zigzag-type AFM state 
in the DFT+I/ {U = 2.75 eV). The obtained spin up (blue) and down (red) densities, respectively, for 2,3,4 and 1,6,5 Ir atoms 
in the hexagon are consistent with the zigzag-type AFM state, (e) The charge density plot of the unoccupied Ir t 2 g band of 
zigzag-type AFM state in the DFT-I-I7 {U = 2.75 eV). (f) The spin density plot of zigzag-type AFM state in the DFT-|-t/ 
(U = 2.75 eV) (different view from (d)). Isosurfaces are chosen as 0.020 (e/a.u.®) for charge densities, and 0.015 (e/a.u.®) for 
spin densities. 


the DFT-I-DMFT calculation was done for the PM state 
(not AFM state) at T = 300 K well above Tj^=lb 
K. Therefore, the gap opens not due to the AFM cell¬ 
doubling but due to the Coulomb correlation effect. This 
result demonstrates that Na 2 lr 03 can be classified into 
a Mott-type insulator. This Mott insulating state of 
Na 2 lr 03 well above Tm is different from the suggested 
marginal Mott insulating state for Sr 2 lr 04 . 1^ It is also 
noteworthy in Fig. [ijb) that, due to the correlation- 
induced imaginary part of the dynamical self-energy, 
t 2 g-driven PDOSs in the DFT-t-DMFT are significantly 
broadened with respect to those in the DFT. 


The correlation-induced band gap and incoherent fea¬ 
tures are more clearly seen in the spectral function plot 
in Fig. (He). As shown in the Supplemental Material, [s^ 
the spectral function becomes more and more incoherent 
with increasing [/, which clearly indicates that the inco¬ 
herent feature is induced by the correlation effect. It is 
also notable that the incoherent feature is more promi¬ 
nent for je//=l/2 than for je//=3/2 states. 


The correlation-induced incoherence in the 
DFT-I-DMFT is essential to describe the measured 
PES of Na 2 lr 03 . In Fig. [2 Ir 5d dominant PES 
near Ep is compared with Ir 5d PDOS obtained from 
both the DFT-bC/ and the DET-bDMET. In the PES, 
there seem to be two dominating peaks, which are 
significantly broadened. 111 Our DET-f-DMET result 


shows the pronounced two peaks that are broadened by 
the correlation-induced incoherence. This dynamical 
correlation-induced incoherence feature is consistent 
with previous report on Na 2 lr 03 ,[ 2 ^ which analyzed 
the PES spectra based on the model Hamiltonian. It 
is also shown in Eig. [2 that, the DFT-|-t/ with static 
correlation yields four peaks, which is not consistent 
with the measured PES. 

When describing the magnetic interaction in the Mott 
insulating state of Na 2 lr 03 , the charge distribution of 
t 2 g-driven unoccupied state is usually considered as a 
rigid object, which is determined by local interactions 
such as the crystal field and the SOC.jl, [3, [3| We 
have confirmed that onset of the zigzag-type AFM or¬ 
dering induces the significant redistribution of the un¬ 
occupied Ir t 2 g state, while the local dynamical cor¬ 
relation in the PM state does not induce the redistri¬ 
bution. Figure [2 illustrates this argument. In Fig. 
[2(b), the charge density of the unoccupied Ir t 2 g band 
in the DFT is plotted. The cubic-shaped charge den¬ 
sity indicates the dominant Je//=l/2 character of the 
band. Figure H c) presents the projection Z{uj,U) be¬ 
tween numerical je//=l/2 bases in the DFT and the 
DFT-I-DMFT with varying U. jeff=i/2 bases were ob¬ 
tained from the diagonalization of the hybridization func¬ 
tion at zero and infinite frequency (A(0) and A(oo)), 
with and without considering the self-energy contribu- 
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FIG. 4: (Color online) Experimental optical condnctivities 
( 7 ( 0 ;)’s of Na 2 lr 03 at T = 300 K, Exp 1 from Ref. [ll[ and 
Exp 2 from Ref. , are compared with calculated optical con¬ 
ductivities using the DFT-I-DMFT, the DFT-|-t/_(?7 = 2.75 
eV), and four-site cluster multiplet calculations. 


tion for the DFT-I-DMFT and the DFT, respectively. 
je//=l/2 orbital for a; = 00 is relevant to the single ion 
anisotropy, while je//=l/2 orbital for a; = 0 is relevant 
to the low energy electronic excitation. 15, SSj As shown 
in FiglHlJc), Z{uj,U) is nearly constant with respect to 
the variation of U, and close to 1 for a; = 00 and 0.995 
for w = 0. This result implies that numerical jeff=il‘2 
bases of the DFT and the DFT-I-DMFT resemble each 
other. Namely, the local correlation does not affect the 
charge distribution, and so the DFT-fDMFT has similar 
charge distribution of unoccupied t 2 g-driven band to that 
of the DFT. 

On the other hand, the onset of the zigzag-type AFM 
ordering induces significant elongation of charge density 
along the ferromagnetic chain direction, as shown in Fig. 
ISKe) for the DFT-|-f/ result. Spin densities in Fig. |3Kd) 
and (f) also verify this phenomenon. In the Je//=l/2 
basis, the spin densities from three different orbitals in 
the t 2 g manifold have equal up, up, and down spin con¬ 
tributions. In our DFT-I-C/ results with the zigzag-type 
AFM ordering, the orbital component directed normal 
to the ferromagnetic chain is much smaller than other 
two orbital components directed along the ferromagnetic 
chain. As a result, the ratio between orbital and spin mo¬ 
ment (/ri//is=1.59 with ^tot=0-70 ^lb) deviates largely 
from the ideal Jeff=l/2 case of ^l/^J'S=‘2■ This feature 
indicates that there will be significant redistribution of 
charge in the AFM transition due to the extended na¬ 
ture of Ir 5d orbital, [ 3 ^ which is not expected for the 
superexchange-based magnetic interaction in the normal 
Mott insulating state. 

In Fig. m experimental optical conductivities tT(w)’s 
are compared with computed optical conductivities from 
the DFT-I-DMFT, the DFT-|-t/, and four-site cluster 
multiplet calculations. ( 2 ^ In experiments, only one peak 
is dominant in the optical spectrum. [ll|, [l7| In the 
DFT-|-[/, there occur two main peaks, as in previous 


DFT-I-I/ calculation,!^ which does not seem to be con¬ 
sistent with experiments. The peak near 0.75 eV should 
be suppressed. In the DFT-I-DMFT, the peak near 0.75 
eV is significantly suppressed due to the incoherency from 
the local self-energy. However, overall peaks are shifted 
upward and the main peak position (^2 eV) is not con¬ 
sistent with the experimental main peak position (~1.5 
eV). If we use smaller U in the DFT-I-DMFT, the main 
peak position is matched but the peak near 0.75 eV is not 
suppressed, and so the resulting a(uj) becomes similar to 
that of the DFT-l-17 (see the SM).[3^ The inconsistent 
DFT-I-DMFT result suggests that the local dynamical 
correlation is not sufficient to describe optical conduc¬ 
tivity of Na 2 lr 03 , which corresponds to a two-particle 
property. 

In contrast, the optical conductivity from the four- 
site cluster multiplet calculation seems to be well con¬ 
sistent with experiments, suggesting the importance of 
the non-local correlation effect in Na 2 lr 03 .[ 2 ^ Simi¬ 
lar feature was found in cuprate too, for which PES 
that is single particle property is well described by the 
single-site DMFT, but optical conductivity that is two- 
particle property is described only by the cluster DMFT 
calculation. ^ Thus, it is expected that the non-local 
correlation effect would be important for optical con¬ 
ductivity due to the extended nature of Ir 5d orbital in 
Na 2 lr 03 . More systematic study of cluster-size depen¬ 
dency of the optical spectrum would be an interesting 
future subject. 

In conclusion, we show that Na 2 lr 03 is a Mott-type 
insulator. Differently from Sr 2 lr 04 , the Mott insulating 
state persists in the PM state well above T^. The local 
correlation-induced incoherence in this Mott insulating 
phase explains the measured PES spectrum of Na 2 lr 03 
well. Yet, due to the extended nature of 5d orbital, the 
insulating nature of Na 2 lr 03 is different from that of a 
standard Mott insulator. The onset of AFM ordering in¬ 
duces significant redistribution of charge. The analysis 
of optical conductivity suggests that the non-local cor¬ 
relation effect also plays a role in Na 2 lr 03 . Therefore, 
our results indicate that Na 2 lr 03 has the Mott-type in¬ 
sulating nature, but the itineracy of the charge degree of 
freedom and the non-local correlation effect should also 
be taken into account to describe physical properties of 
this system having localized and itinerant duality of 5d 
electrons. 
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